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Abstract 

We report the first observation of the decay modes B° — ► D + D~ , B~ — ► D°D~ and B~ — ► 
L>°L>*- based on 152 x 10 6 55 events collected at KEKB. The branching fractions of 5° -> D + D~, 
B~ -► L>°L>- and B~ -► D°L>*- are found to be (3.21±0.57±0.48) x 10" 4 , (5.62±0.82±0.65) x 
and (4.59 ± 0.72 ± 0.56) x 10 -4 , respectively. Charge asymmetries in the B~ — ► D°D~ and 
i? _ — ► D°D*~ channels are consistent with zero. 

PACS numbers: 13.25.Hw, 11.30 Er 
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Mixing induced CP-violating asymmetries in b — > ccs decays have been observed at the 
B factory experiments, Belle and BaBar, at levels consistent with Standard Model (SM) 
predictions Cabibbo-suppressed double charm decays (e.g. B° — > _dM+_D(*) _ ) are 
dominated by b — ► ccd tree diagram contributions. Additional penguin contributions with a 
different weak phase are expected to be small in comparison with tree diagram contributions. 
Hence, time dependent CP-violating asymmetries in such double charm decays should be 
nearly equal to those in B° — ► J/ipKs {b — ► ccs) type decays. However, a variety of 
processes beyond the SM can provide additional sources of CP violation |2|. Thus, the 
B° _> £>(*)+ £>(*)- decay modes can be used to confirm the SM predictions of CP violation 
or to look for physics beyond the Standard Model. 

Signals for the decay modes B° -> D* + D*-, B° -> and P° -> £)+£)*- have 

already been observed. Measurements of time dependent CP-violating asymmetry parame- 
ters are also available for these modes [il]. While the branching fraction for B° — > D + D~ is 
expected to be fairly large, evidence for this decay mode has not yet been reported. Using 
SU(3) symmetry and the world-average branching fraction jfjj], B(B° — ► D + D~) is estimated 
to be sin 2 9 C x B(B° -> £>+D~) ~ (4.0 ± 1.5) x 10~ 4 , where fl c is the Cabibbo angle. 

The decay modes B~ — > D°D~ and P~ — > D°D*~ are expected to be dominated by 
tree diagrams with some additional contributions from penguin diagrams. Scaling from 
the well-measured branching fractions for the Cabibbo-favored processes B~ — > D°D~ and 
B~ — > D°D*~/B~ — > D*°D~, one can estimate the branching fractions of these two decay 
modes to be (6.5±2.0) x 10 -4 and (5.1±2.2) x 10~ 4 , respectively. Assuming SU(3) symmetry, 
measurement of these branching fractions will enable better understanding of the penguin 
processes. 

Here, we report the first observation of the decays B° — ► D + D~, B~ — > D°D~ and 
B~ — > D°D*~ . Inclusion of charge conjugate modes is implied throughout this paper. The 
analysis is based on a 140 fb" 1 data sample at the T(4S) resonance (10.58 GeV) and a 
16 fb -1 data sample 60 MeV below the T(4S) peak (referred to as off-resonance data), 
collected with the Belle detector at the energy asymmetric e + e~ collider KEKB 0. The 
data sample contains 152 xlO 6 BB events. The fractions of neutral and charged B mesons 
produced in T(4S') peak are assumed to be equal. 

The Belle detector is a general purpose magnetic spectrometer with a 1.5 T magnetic field 
provided by a superconducting solenoid. Charged particles are measured using a 50 layer 
Central Drift Chamber (CDC) and a three layer double sided Silicon Vertex Detector (SVD). 
Photons are detected in an electromagnetic calorimeter (ECL) consisting of 8736 CsI(Tl) 
crystals. Exploiting the information acquired from an array of 128 time-of-flight counters 
(TOF), an array of 1188 silica aerogel Cerenkov threshold counters (ACC) and dE/dx- 
measurements in the CDC, we derive particle identification (PID) likelihoods C^IK- A kaon 
candidate is identified by a requirement on the likelihood ratio Lk/{£k + £vr) such that the 
average kaon identification efficiency is ~ 93% with pion misidentification rate of ~ 10%. 
Similarly, charged pions are selected with an efficiency of ~ 95% and kaon misidentification 
rate of ~ 10%. We select charged pions and kaons that originate from the region dr < 1 cm 
and \dz\ < 4cm with respect to the run dependent interaction points (IP), where dr, dz are 
the distances of closest approach of ir/K tracks to the IP in the plane perpendicular to and 
along the z-axis (the z-axis is defined as passing through the nominal interaction point and 
antiparallel to the positron beam). All tracks compatible with the electron hypothesis (~ 
0.2% misidentification rates from pion/kaon) are eliminated. No attempt has been made to 
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identify muons, which represent a background of about 2.7% to the pion tracks. 

Neutral kaons (Ks) are reconstructed via the decay K$ — * 7i + tt~ with no particle iden- 
tification requirement for daughter pions and the two-pion invariant mass is required to be 
within 11 MeV/c 2 (~ 3.5cr, where a is the invariant mass resolution of tt + 7t~) of the K$ 
mass. To improve the purity of K$ candidates, we impose K$ momentum-dependent criteria 
on the impact parameter of the pions, the distance between the closest approaches of the 
pions along the beam direction, the distance of the 7i + n~ vertex from the interaction point, 
and the azimuthal angle difference between the direction of 7r + 7r~ vertex from the IP and 
the Ks momentum direction. Mass and vertex constrained fits are applied to obtain the 4- 
momenta of Ks candidates. Neutral pions (tt°) are reconstructed from pairs of isolated ECL 
clusters (photons) with invariant mass in the window 118 MeV/c 2 < M 77 < 150 MeV/c 2 
(~ ±3cr). The energy of each photon is required to be greater than 30 MeV in the barrel 
region, defined as 32° < 6* 7 < 128°, and greater than 50 MeV in the endcap regions, defined 
as 17° < # 7 < 32° or 128° < 9 7 < 150°, where 6 1 denotes the polar angle of the photon. 
Mass constrained fits are applied to obtain the 4-momenta of 7r° candidates. 

Beam gas events are rejected using the requirements \P Z \ < 2 GeV/c and 0.5 < E vis / \J~s < 
1.25, in the T(4S) rest frame, where P z and E vis are the sum of the longitudinal momentum 
and the energy of all reconstructed particles, respectively, and \fs is the sum of the beam 
energies in the T(4S) rest frame. The continuum (e + e~ — > qq, where q = u,d, s, c) events are 
suppressed by requirements on the ratio of the second to the zeroth Fox- Wolfram moments 
@, R 2 < 0.35. 

The D° meson is reconstructed through its decay to 
K~n + , K~tt + 7t°, K~7t + 7t + tt~ , K sit + -k~ and K + K~ . The D + meson is reconstructed 
through its decay to K~7t + tt + and K~K + tc + . Mass and vertex constrained fits are applied 
to improve the D meson momentum resolution. 

The tracks from the D decays are chosen with a criterion, \dz^ — dz^\ < 2 cm, where dzm) is 
the distance of closest approach of track to the IP along the z-axis. Large combinatorial 
backgrounds are removed by requiring that the invariant mass of daughter particles is within 
2.5 a from the nominal Z)-mass, where a (~ 4.5 MeV/c 2 ), the mass resolution, depends on 
the decay chain. 

D* + candidates are reconstructed by combining the D° with a slow charged pion with 
dr < 2 cm and \dz\ < 10 cm with respect to the D vertex. D° mass windows are widened 
to ± 20 MeV/c 2 for the reconstruction of D* + candidates. D* + candidates are required to 
have a reconstructed mass difference between the D* + and D° within 2.0 MeV/c 2 (~ 3.0cr) 
of the nominal mass difference. A kinematic fit with the D* + mass is applied to obtain the 
4-momenta of the D* + candidate. 

To reduce large combinatorial backgrounds in D° decays to K~7r + 7i + n~ , a tighter impact 
parameter requirement, dr < 0.5 cm, is applied to all four tracks. Note that this is not ap- 
plied for D* + — > D°(—>- K~7t + 7t + 7t~)tt + signals because this mode has much less background. 
The D° — > K _ 7T + 7r° decay mode is used only for D* + candidates to avoid large backgrounds 
in D° signals in this decay mode. 

Combinations of DD^ are used to reconstruct candidate B mesons. DD^ signals are 
contaminated with background from misidentified D mesons and also combination of two D 
candidates from opposite B mesons. There are two important kinematic variables to extract 
signals from these backgrounds, (i) the energy difference, AE, between the measured energy 
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of the candidate event and the beam energy, -Ebeam, m the T(4S f ) rest frame and (ii) the 



beam energy constrained mass, M^ c = y E^ clim — (£\ Pi) 2 , where are momentum vectors 

of the primary £)W candidates. The AE distribution is used to extract the signal yield 
since peaking backgrounds are expected in the Mt, c distributions. The fit is performed for 
events where M bc satisfies 5.272 GeV/c 2 < M bc < 5.288 GeV/c 2 and the fit range in AE 
is from —70 MeV to 200 MeV. The restricted range in negative AE is chosen to exclude 
contributions from other B decays, such as B° — > D°D* (D*° — > D°j). Selected events 
contain multiple B candidates with a multiplicity depending on the signal channels, which 
varies from 1.02 to 1.07. In events with more than one candidate B meson, the candidate 
with the smallest x 2 (= {AM Dl / a Mr) ^) 2 + (AM D2 /<j Md ^) 2 ) is chosen, where AM D . is the 
difference of the reconstructed and nominal mass of Dj, and cr A/o is the resolution in Mo r 
In B~ — > D°D*~ candidates, a (AMp- / Va* _) 2 term is also added to the x 2 to choose 
the best candidate, where AMd*- is the difference between measured and nominal mass 
difference between D*~ and D° , and ctm d ,_ is the resolution in measured mass difference. 

Unbinned extended maximum likelihood fits to AE distributions are used to extract the 
signal yields. The signal shape is modeled as a sum of two Gaussians, 

/ \ I- — n \~ 
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which is dominated by a Gaussian of width, a ~ 6 MeV and a wider Gaussian, whose width 
is ~ 2-3 times larger than the main Gaussian function, but whose contribution is only 
15%— 17% of the main Gaussian function. Backgrounds are modeled with a linear function. 

Signal Monte Carlo (MC) is fitted with this function to determine f\ and The fitted 
values of fx and fi are 0.051 and 2.48, respectively. The AE distribution in data is wider 
than in MC. The decay channel, B~ — > D°D~ is used as a control sample to calculate the 
scaling factor, / ( = [<Jdata/ ^mc] b-->d°d~ = ^ 0-05 ) of the Gaussian width in data. 

To extract signal yields from data, the width of the main Gaussian function is fixed to 
v — /cmc an d remaining four parameters (parameters A and \x of signal function and two 
parameters of the linear background function) are determined in the fit. 

Signal and backgrounds are studied with Monte Carlo event samples that are generated 
using the QQ event generator The response of the Belle detector is simulated by a 
GEANT3-based program ^(J. The simulated events are reconstructed and analysed with 
the same procedure as is used for the real data. A large generic BB MC samples, with 
a luminosity equivalent to 330 fb _1 of data is used to look for peaking backgrounds in 
AE distributions. We have also studied feed across among these signals and from other 
B — > D^D^ decays. Background due to continuum events are studied by analysing the 16 
fbT 1 of off-resonance data and simulated MC events equivalent to ~ 330 ftT 1 of data. No 
peaking backgrounds are observed in AE distributions of these samples. Signal efficiencies 
in B° -> D+D-, B~ -> D°D~ and B~ -> D°D*- decay modes are 11.9±0.1%, 11.1±0.1% 
and 4.8±0.1%, respectively. 

Figure ^ shows the AE distributions in data for the decay modes B° — > D + D~, B~ — > 
D°D~ and B~ — > D°D*~ . There are clear structures near AE = 0. The results of the fits for 
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FIG. 1: Fit results of AE distributions in data. Points with error bars are the observed events in 
data, solid lines are the results from the fit and dashed lines represent the background components. 



TABLE I: Observed signal yields, statistical significances (a) and branching fractions. 



Channel 




iVobs 




a 






B x 10 4 






B° -» D + D- 


54.3 


± 


9.7 


7.3 


3.21 


± 


0.57 


± 


0.48 


B~ D°D- 


120.5 


± 


17.6 


8.0 


5.62 


± 


0.82 


± 


0.65 


B~ -» D°D*- 


73.6 


± 


11.5 


8.2 


4.59 


± 


0.72 


± 


0.56 



these decay modes are also shown in these plots. The signal yields obtained from the fits are 
given in Table UJ The statistical significance of the yields, defined as \f— 2 In (£ /£ maiX ), is 
7.3, 8.0 and 8.2 for the B° -> D+D~, B~ -> D°D~ and B~ -> D°D*~ channels, respectively, 
where £ (£ max ) is the maximum likelihood without (with) the signal contribution. The 
corresponding AE distributions for events in the M^ c sideband region (5.2 GeV/c 2 < M^ c < 
5.26 GeV/c 2 ) are also checked and they do not show any structure. To check for possible 
background from modes such as B — > DKtt, B — > DKittt or charmless 5 — > KirK(n7r), 
we also examine the 5 signal yield for combinations when one of the D candidates has an 
invariant mass in a sideband outside the nominal D mass window. No significant yield is 
observed in such combinations. 

Branching fractions obtained for these three modes are listed in Table HI where the first 
error is statistical and the second error is systematic. This is the first measurement of the 
branching fractions for these decay modes. All results are consistent with the expectation 
from SU(3) symmetry. As a consistency check, the M^ c distributions are also fitted after 
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constraining \AE\ < 40 MeV and they show consistent signal yields. 

The distribution of the helicity angle for B~ — > D°D*~ channel is also studied. The 
helicity angle, 0, is defined as the angle between the direction opposite to the B meson and 
that of the slow pion in the the D*~ rest frame. Figure 121 shows the AE sideband subtracted 
(the signal region is defined as |Ai?| <20 MeV and sideband regions are 50 MeV < |A.E| < 
70 MeV) helicity angle distributions in data and for the MC signal. The data follow a cos 2 
distribution as expected for a B to pseudoscalar-vector decay. 




-1 -0.5 0.5 1 

COS© 



FIG. 2: Background-subtracted distribution of the cosine of the helicity angle in D*~ —* D ir~ 
decay for B~ —* D°D*~ candidates in data (points with error bars) and for B~ —* D°D*~ signal 
MC (hatched histogram). 

The systematic uncertainty, shown in Table |Hj is obtained from a quadratic sum of the 
uncertainties in (a) the track finding efficiency, ranging from 1% for high momentum tracks 
to 8% for pions of 80 MeV/c, estimated from partially reconstructed D*~ — > D°(— > K$(^ 
7t + 7t~)7t + 7t~)tt~ events and a track embedding study; (b) the ir° reconstruction efficiency, 
estimated from a comparison of D° — > K~7c + tt° yields in data and MC ; (c) the K$ selection 
efficiency, estimated from a comparison of D° — > Kstt + it~ yields in data and MC; (d) K/ir 
selection efficiencies, estimated using D*~ — > D°(— > K + tt~)tt~ events; (e) The world-average 
D°,D + and D* + branching fractions |5j; (f) scaling factor for AE distributions in data; 
(g) MC statistics; (h) the total number of BB events {N BB ). The systematic error is also 
studied by (a) varying the AE fit ranges within —200 MeV to +200 MeV; (b) the choice 
of the fitting functions; (c) deriving the branching fraction without selection of the best 
candidate. Systematic uncertainties from the latter three sources are negligible compared 
to those from the other sources. 

Charge asymmetry, A = (JV_ -N + )/(N_ + N + ) in B~ ->■ D°D~ and B~ -> D°D*- 
channels is —0.05 ± 0.15 ± 0.05 and 0.15 ± 0.15 ± 0.05, respectively, where N_(N + ) is the 
number of observed events in B~(B + ) decays. Systematic errors on charge asymmetries are 
determined from high statistics B~ — ► D and B — > D*(mr) event samples. 

In summary, we report first observations of the decays B° — > D + D~, B~ — > D°D~ and 
B~ — > D°D*~ using 152 million BB events. We measure branching fractions for these three 
decay modes, which are consistent with the expectation from SU(3) symmetry. Charge 
asymmetries in the B~ — > D°D~ and B~ — > D°D*~ are consistent with zero. We find a 
statistically significant signal in the B° — > D + D~ channel; in the future this mode will be 
used to perform additional studies of time dependent CP violation in b — > ccd decays. 

We thank the KEKB group for the excellent operation of the accelerator, the KEK cryo- 
genics group for the efficient operation of the solenoid, and the KEK computer group and the 
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TABLE II: Individual contributions of systematic uncertainties (in %). 





D+D- 


D°D- 


D°D*- 


Track finding 


6.5 


7.2 


9.1 








2.0 


K s 




0.4 


1.0 


K/ir selection 


5.5 


5.5 


5.1 


AE scale 


2.4 


2.6 


2.2 


MC statistics 


0.9 


1.2 


2.0 




0.5 


0.5 


0.5 


D branching fractions 


12.1 


6.6 


5.3 


Total 


15.0 


11.5 


12.1 
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